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ABSTRACT 

Aims. We wished to analyse a sample of observations from the XMM-Newton Science Archive to search for evidence of exospheric 
solar wind charge exchange (SWCX) emission. 

Methods. We analysed 3012 observations up to and including revolution 1773. The method employed extends from that of the previ- 
ously published paper by these authors on this topic. We detect temporal variability in the diffuse X-ray background within a narrow 
low-energy band and contrast this to a continuum. The low-energy band was chosen to represent the key indicators of charge exchange 
emission and the continuum was expected to be free of SWCX. 

Results. Approximately 3.4 % of observations studied are affected. We discuss our results with reference to the XMM-Newton mis- 
sion. We further investigate remarkable cases by considering the state of the solar wind and the orientation of XMM-Newton at the 
time of these observations. We present a method to approximate the expected emission from observations, based on given solar wind 
parameters taken from an upstream solar wind monitor We also compare the incidence of SWCX cases with solar activity. 
Conclusions. We present a comprehensive study of the majority of the suitable and publically available XMM-Newton Science 
Archive to date, with respect to the occurrence of SWCX enhancements. We present our SWCX-affected subset of this dataset. The 
mean exospheric-SWCX flux observed within this SWCX-affected subset was 15.4 keVcm"^ s"' sr"' in the energy band 0.25 to 
2.5 keV. Exospheric SWCX is preferentially detected when XMM-Newton observes through the subsolar region of the Earth's mag- 
netosheath. The model developed to estimate the expected emission returns fluxes within a factor of a few of the observed values in 
the majority of cases, with a mean value at 83 %. 

Key words. X-rays - diffuse background - Sun: solar-terrestrial relation - methods: data analysis 



1. Introduction 

This paper follows that of ICarter & Sembavl (l2008h (hereafter 
Paper I), whereby a set of approximately 180 XMM-Newton ob- 
servations, taken between revolutions 52 and 1 104 (March 2000 
until December 2005), were analysed to search for cases of solar 
wind charge exchange emission (SWCX) occurring within the 

' Earth's magnetosheath or in near interplanetary space. In Paper 
I we searched for time-variable SWCX signatures and found 
that approximately 6.5 % of the observations studied were af- 
fected. One case of SWCX showed an extremely rich emission 
line spectrum, which was attributed to a passing C oronal Mass 

' Ejection (CME) and has been discussed in detail bv ^Carter et al.l 
(12010 ). This paper extends the work of Paper I, to provide a 
complete sample, covering 3012 suitable XMM-Newton obser- 
vations downloaded from the XMM-Newton Science Archive 
(XSAfl We include data from between revolutions 28 and 1773 
(February 2000 until August 2009). 

SWCX processes occur at many locations within the so- 
lar system including planetary magnetosheaths, the corona of 
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comets, within the heliosphere and at the heliospheric bound- 
ary where the outer reaches of the interplanetary magnetic field 
encounter that of the surrounding interstellar medium. Those 
charge exchange processes that result in the emission of X- 
rays occur for interactions involving a highly-charged solar wind 
ion, for example the bare oxygen ion O^^, and a donor neu- 
tral species, such as hydrogen in the case of geocoronal SWCX 
emission. 

The possibilities for viewing exospheric SWCX emission in 
the vicinity of the Earth depend on the orbital and viewing con- 
straints of the observatory in use, however one may be able to 
detect enhancements due to charge exchange occurring within 
the heliosphere. This will show fluctuations on long er timescales 
to SWCX occurring within the Earth's exosphere (ICravens et al.l 
2001). SWCX emission from within the heliosheath (for exam- 
ple resulting from the helium-focusing cone) is a contributor 
to the X-ray emission from the supposed Local Hot Bubble in 
which the Sun resides (Ko utroumpa et al] 2007, 2008), and to 
what level it contributes is under great debate. We concentrate, 
however, on SWCX emission occurring in the near vicinity of 
the Earth and use the time variable nature of this emission as our 
marker for selecting affected XMM-Newton observations. 

In addition to the pioneering work noting the so-called 
Long Term Enhancements (subsequently know n to be due to 
SWCX emission) within ROSAT observations dSnowden et al.l 
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Il995h . the assignment of such enhancements to SWCX emission 
has also been observed in studies undertaken using data from 
Suzaku (Fujimoto et al. 2007; Bautz et al. 2009) and Chandra 
dSmith etalJl2005i) . Enhancements in soft X-ray band XMM- 
Newton spectra have been attributed to SWCX contamination 
in the literature jSn owden et al. 200 4 iKuntz & Snowdenll2008t 
ISnowdenetani2009t iHenlev & S helton"2008'). This has almost 
exclusively involved the comparison of multiple pointings of the 
same field which has enabled the serendipitous detection of the 
low-energy enhancement, most notably around the O vii helium- 
like triplet at approximately 0.56 keV. Observations of the Groth- 
Westfall Strip, Polaris Flare region and the Hubble Deep Field 
as discussed in Kuntz & Sno wden (2008) were included in the 
analysis of Paper I as control cases f or our method. 

The XMM-A^evvfon observatory (iJansen eta"ni2001h has the 
largest collecting area in the band 0.2 to lOkeV of all X- 
ray telescopes currently in orbit and the European Photon 
Imagi ng Camera (EPIC) suite of in struments on board (two 
MOS (iTurner et al.ll200ll) and one pn (IStriider et al.ll200lh CCD 
cameras) provide medium spectral resolution ~ 17). During 
various periods of its orbit and depending on observational con- 
straints, XMM-Newton may view regions of the Earth's magne- 
tosheath which are predicted to exhibit the highest X-ray emis- 
sivity due to SWCX between highly charged solar wind ions and 
hydrogen in the Earth's exosphere (Robertson et al. 2006, and 
references therein). We use data from the EPIC-MOS cameras 
to identify times of variability in a low energy band that is not 
mirrored in a higher energy band, to select those periods we sus- 
pect have been affected by variable SWCX. 

SWCX emission can be considered as a contaminant 
by sections of the community wishing to study Galactic or 
Extragalactic sources beyond the heliosheath. Understanding the 
occurrence and nature of the emission is important therefore 
for the process of elimination from astronomical observations. 
However, SWCX emission can be used to provide diagnostics 
of the solar wind. X-rays emitted from the comas of comets 
have been suggested and used for this purpose JCra vens 1997; 
iDennerletalTI 19971: iLisse et alj 120011: iBodewits et al.i i2007D . 
A detailed study of the heavy ion constituents of a CME 
was discussed in ICarter et"an (l2010l) . In addition, exospheric- 
SWCX may be used to test models of dynamical processes 
within the Earth's magnetosheath and the solar-terrestrial con- 
nection, such as flux-transfer events or boundary layer phenom- 
ena (Collier et al. 2010). 

The layout of the paper is as follows. In Section |2] we de- 
scribe the method employed to identify those observations of 
interest. In Section[3]we present the overall results for the whole 
sample studied and compare the occurrence of SWCX with the 
solar cycle. In Section |4] we describe fitting spectral models to 
cases of SWCX enhancement data. In Section |5] we discuss in 
more detail several of those observations affected by SWCX. In 
Section |6] we describe a model of the expected X-ray emission 
using the orbit of XMM-Newton and the solar wind conditions 
at the time of the observation. We finish with our discussion and 
conclusions in Section|2l 

2. Data analysis 

All observational data used in this study are available through the 
XSA where we extract the Original Data File s (ODF). We con- 
sider observations taken by the EPIC-MOS (ITurner et alJl200lh 
cameras (MOSl and MOS2) using the full-frame mode only. 
Therefore the observations used provided an even sample across 
the mission and the event list data from the cameras could be 



combined if the same filter was used for both. We considered 
observations up to and including revolution 1773 (August 2009). 

Individual observation event lists were created from the 
Original Data Files for each EPIC-MOS camera and filtered for 
soft-proton contamination using the publicly available Extended 
Source Analysis Software (ESASj^ Package (mos-filter tool). At 
the time of data processing, ESAS was only available for the 
EPIC-MOS cameras. This tool fits a Gaussian to a histogram of 
in-field-of-view count rates in a high energy band (2.5 keV to 
12.0 keV). Time periods with count rates beyond a threshold of 
+ 1.5 cr away from the peak of this Gaussian were removed by 
applying a Good Time Interval (GTI) file to the event lists. The 
GTI files created for each EPIC-MOS event list were combined 
together to form one EPIC-MOS GTI file and this file was reap- 
plied to both the MOS 1 and MOS2 event lists to provide simulta- 
neous coverage during each observation. Re solved point sources 
(from lists created for the 2XMM catalogue (IWatson et al.ll2009h 
using a minimum likelihood threshold > 6), were removed from 
the field of view by extracting events in a circular region of 35 
arcseconds radius about the source position. Those observations 
judged, after a visual inspection (and prior to creation of any 
spectra, see Section |4]i, to show residual source contamination 
(the wings of the point spread function of the EPIC-MOS cam- 
era being evident in an image of the event file) passed through 
an additional spatial filtering stage using a larger extraction re- 
gion to further clean the dataset. A more detailed description of 
the filtering steps an d nature of soft prot on contamination can be 
found in Paper I and lCarter et al.l (l2010h . 

The method employed here followed the first steps as de- 
scribed in detail within Paper I. In summary, two lightcurves 
with bin size of 1 ks were created for each observation from 
events within the full field-of-view (radius of 13.3 arcminutes). 
When both EPIC-MOS cameras were used during an observa- 
tion and employed the same filter, events from both cameras 
were used to construct the lightcurves. The first lightcurve was 
chosen to represent the continuum, covering events with energies 
in the range 2.5 keV to 5.0 keV. The second lightcurve was ex- 
tracted using events with energies in the range 0.5 to 0.7 keV to 
cover the strong SWCX emission from O vii and O viii (the line- 
band lightcurve). This energy range incorporates emission ener- 
gies from the O vii triplet and resonance lines which are dom- 
inated by the forbidden line transition at 0.56 keV. Lightcurves 
were then exposure-corrected for periods removed during the fil- 
tering steps. As the MOSl and MOS2 event files for each obser- 
vation, prior to the lightcurve creation, were filtered using a sin- 
gle GTI file which is not energy specific, the exposure-coverage 
for each bin was the same for the line-band and continuum 
lightcurves. Therefore the same exposure-correction factor for 
an individual bin was applied to both lightcurves in this step. We 
keep bins of the lightcurve with at least 60% of the full exposure 
for that bin and reject the remaining bins (in contrast to Paper I 
when we kept all bins with at least 40% exposure). Lightcurves 
were rejected from further analysis if they were less than 5 ks in 
length. By increasing the strictness of the bin coverage thresh- 
olds, we reduced the incidence of type I errors (those incorrectly 
labelled detections), but ran the risk of increasing the number 
of cases that were incorrectly labelled as non-detections (i.e. not 
showing a deviance from the null hypothesis of a linear fit be- 
tween the line-band and continuum (type II errors). We scaled 
each lightcurve by its mean to produce adjusted lightcurves. The 
count rates for the line-band and the continuum band were then 
always of the same order which facilitated the identification of 
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periods of enhancement in the line-band lightcurve. An exam- 
ple of the combined-MOS lightcurve adjustment process can be 
found in Figure [1] (panels top-left, top-right and bottom-left). 

We plotted a scatter plot between the two bands (using the 
adjusted line-band as the dependent variable), shown in Figure[T] 
(bottom-right). A linear model fit to each scatter plot was com- 
puted using the IDL procedure. Unfit, which minimises the 
statistic. 

We computed the reduced-;^'^ for the fit, hereafter referred to 
as by dividing the x^ by the number of bins minus one, to 
account for the reduction in the number of degrees of freedom 
made by fitting a linear model to the data. A high x^ indicates 
that a significant fraction of the points deviate significantly from 
the best fit line. We expected these cases would be more likely 
to show variable SWCX-enhancement. In addition we computed 
the x^ values for each individual lightcurve in terms of the de- 
viation from the mean of that lightcurve. We calculated the ratio 
between the line-band and continuum x^ values to add to our 
diagnostic (hereafter denoted as /?^). 

3012 observations made up the final sample to be used 
for further analysis. The results from Paper I showed us that 
those observations exhibiting both high x]j and high were 
most likely to show near-Earth time-variable SWCX signatures. 
Observations that fulfilled these criteria were considered for fur- 
ther analysis, spectrally, temporally and with regard to the orien- 
tation of XMM-Newton. 

3. Global results 

All observations in our final sample, after any rejections as de- 
scribed below, were ranked by x^- The two highest ranked ob- 
servations were observations of comets. Although resolved point 
sources have been removed, cometary X-rays are diff'use and will 
likely be spread over a large fraction, if not all of the field of 
view. We assume that the dominant variations in the line-band 
lightcurve that result in such high ranking are due to SWCX 
emission occurring within the cometary coma and not to any 
emission occurring within the vicinity of the Earth. We discuss 
the cometary cases in Section 1331 

After ranking the observations, we were able to study those 
that exhibited the highest xj, and in more detail on a case by 
case basis. 

Observations were examined for residual point or extended 
sources that may contribute to the high variations in the line- 
band. Extended or diff'use residual sources that remain in the 
field of view will not aff'ect this detection method providing no 
inherent variation occurs within these sources, as expected for 
sources outside the solar system, in either one of the bands. 
Cases were also examined for residual soft proton contamina- 
tion. Although the files used in this analysis have been filtered 
for periods of soft proton flaring, residual contamination may 
remain. Excessive scatter due to residual soft proton contamina- 
tion impedes the ability to identify periods of exospheric SWCX 
and will result in some type II errors in our sample. Also, if exo- 
spheric SWCX occurs throughout the entirety of the observation 
little or no variation can be seen in the line band, resulting in a 
rejection of this case. Excessive and simultaneous variations in 
both the line-band and continuum will result in a high x'^ yet a 
low value of 7?^. We concentrate our analysis therefore on cases 
that exhibit both high;if^ and high 7?^. 

For some observations with short lightcurves, it was impos- 
sible to identify any time-periods of boosted line-band emis- 
sion (the putative SWCX enhancement periods). These observa- 
tions were disregarded as SWCX-enhancement cases. Also, the 



Table 1. Highest ranked observation by x^- Those observations 
with severe residual soft proton contamination have been ex- 
cluded. 



Revn. 


Obsn. 


^« 


Comment 


0369 


0103461 lOI 


856.7 


Comet C2000 WMl (LINEAR) 


0808 


0I64960I0I 


226.6 


Comet C2001 Q4 (Neat) 


0342 


0085150301 


27.2 


Paoer I & Carter et al. (2010") 


0209 


0093552701 


23.0 


Paper I 


1014 


0305920601 


15.0 


Paper I 


0690 


0I4963030I 


14.1 


Paper I 


0623 


0I506I0I0I 


13.5 


New case 


1177 


0406950201 


13.3 


Comet 73p 


0339 


0054540501 


13.2 


New case 


0422 


01 13050401 


12.7 


New case 



method described in this paper tested only for variable SWCX 
on short timescales. Spectral analysis of suspected SWCX cases 
was therefore only possible when a clear line-band enhancement 
period could be identified during the duration of the observation. 

We find 103 observations in our sample that show indica- 
tions of a time-variable exospheric SWCX enhancement, that 
are not excluded from consideration based on the reasons de- 
scribed above. All of these cases had a x^ value greater than 
or equal to 1.2 and a 7?^ value of greater than or equal to 1.0. 
These cases make up only ~20% of all observations that have 
values of x^ and above these thresholds, indicating that al- 
though the values of x]^ and 7?^ are indicators of a SWCX- 
enhancement, considerable inspection of an observation on a 
case-by-case basis is still required. The majority of cases in the 
whole sample have a R^ value less than 1, indicating that there 
is more variation seen in the continuum compared to the line- 
band. The continuum incorporates the break energy at ~3.2keV 
in the two-power law models of residual soft proton contamina- 
tion dKuntz & Snowden 2008). For higher intensity soft-proton 
flares, the slope of the power-law becomes flatter Therefore the 
higher energy and wider continuum will have a greater variance 
than the softer, narrower line-band due to the presence of unfil- 
tered residual soft protons. 

The top 10 observations as ranked hy x^ are given in Table [T] 
(excluding those observations that had been rejected from con- 
sideration). Three of these top ten result from cometary obser- 
vations, four were identified in the work presented in Paper I 
and three are new observations discovered during this analysis. 
Many of the highest ranked cases had previously been identified 
in the literature. It is our intention to make the full ranked list 
of observations used in this analysis publically available online 
in the future, most probably through the XMM-Newton EPIC 
Background Working Group (BGWG)El We include a summary 
table listing the SWCX cases in Table lA. II 

A scatter plot ofx^ versus 7?^ is given in Figure|2] Those ob- 
servations exhibiting both high x]j and high 7?^ but which do not 
show clear SWCX signatures (i.e. an enhancement in the mean- 
adjusted line-band lightcurve compared to that of the continuum) 
have been rejected from further analysis. We discuss the detec- 
tion methods of SWCX cases in the literature, other than by a 
search for a time-variable low-energy component in Section [J!2] 
Cometary X-ray emission is discussed in Section [33] 

In Figure [3] we plot the total number of observations and the 
fraction of observations that show SWCX enhancements, ver- 
sus the GSE-X position of XMM-Newton at the mid-point of 
each observation. We can see from this figure that XMM-A^ewfo« 

^ http://xmm.vilspa.esa.es/external/xmm_sw_cal/background 
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Fig. 1. Lightcurve correction procedure example for observation with identifier 0150680101 (line-band (black), continuum (red) for 
panels top-left and bottom-left). Top left: example lightcurves showing a peak in the line-band that is not reflected in the continuum. 
Top right: exposure coverage for each bin, the threshold at 60% is marked by the red dashed line. Bottom left: lightcurves after the 
adjustments for exposure correction and scaling by the mean. Bottom right: example scatter plot for this observation. 




200 400 600 800 
No. obsn 



Fig. 2. xj, versus R^^. Red filled circles indicate those observa- 
tions with time-variable SWCX signatures. Histograms of the 
xj, and Rj^ values are shown in the side panels (using a bin size 
of 0.25). 



is preferentially found on the subsolar side of the Earth when 
SWCX enhancements occur In addition, the largest bin of the 
fractional plot occurs around the nominal magnetopause stand 



off distance, at approximately lOR^. The XMM-Newton line- 
of-sight for these cases traversed the subsolar region (sunward 
side) of the magn etosheath, as expec ted according to the mod- 
elling work of Rob ertson et al.l ( l2006l) . XMM-Newton viewing is 
constrained by the fixed solar panels and limits imposed to avoid 
directly observing the Sun, Moon and Earth and is only able to 
view the sunward side of the ma gnetosheath at certain times of 
the year (ICarter & Sembavll2008l) . 

We also consider the seasonal variation of the occurrence 
of SWCX in this sample. 64 SWCX cases occurred during the 
summer months (April until September inclusive) and 39 during 
the winter (October until March). More exospheric SWCX cases 
are expected in the summer months, as viewed by XMM-A^ewfo« 
(■Carter & Sembav 2008) . 

3. 1. Relationship witti fhe solar cycle and solar wind 

In Figure|4]we plot the solar sunspot numbeiQ from the latter half 
of Solar Cycle 23 and mark the times at which XMM-Newton 
observations with known SWCX (this paper, paper I and XMM- 
Newton exospheric-SWCX cases in the literature) occurred. We 
plot a histogram of the number of SWCX cases per half year, 
to remove any bias resulting from the seasonal constraints on 
pointing angle experienced by XMM-A^ewfon. Each histogram 
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Fig. 3. Total number of observations (black) versus GSE-X po- 
sition and the fraction of observations detected with exospheric 
SWCX enhancements (red). 




2010 



2004 2006 2008 2010 
Year 

Fig. 4. Top panel: sunspot number versus time. Bottom panel: 
the coloured histogram of the fraction of observations affected 
by exospheric SWCX is binned into six month periods (blue - 
summer, black - winter). The total number of all observations 
for each period is noted above the bin. 



bin starts at the start of summer (!*' April) or start of winter 
(1" October). As expected for cases of exospheric SWCX, there 
are more cases at times of high solar activity around solar maxi- 
mum than when approaching solar minimum. Also, for the cases 
from 2002 onwards and approaching solar minimum, we see a 
higher proportion of SWCX cases in the summer six-month pe- 
riod compared to the winter period, for the same year. The sig- 
nificance of this trend however should not be overstated due to 
the low number of cases in each bin. 



3.2. Multiple pointings of target fields 

Multiple pointings towards the same target allow one to compare 
diffuse and extended emission spectral models that may exhibit 
spectral variations indicative of SWCX contamination. The long 
term enhancements of the ROS AT all-sky maps, which were sub- 
sequently attributed to SWCX, were first identified by compar- 
isons between fields (Snowden et al. 1995). Kuntz & Snowden 
( |2008|) examined multiple XMM-Newton observations of the 
Hubble Deep Field, amongst other targets, and identified a pro- 



porti on of their set affec ted to different degrees by SWCX emis- 
sion. Bautzet an (l2009h inferred SWCX-enhancements in obser- 
vations by Suzaku towards the cluster Abell 1795 after examina- 
tion of a low-energy lightcurve revealed peaks coincident in time 
with enhance ments in the solar w ind proton flux as measured 
by ACE. He nley & SheltonI (1201 Ol) used a large set of XMM- 
Newton observations to compare intensities from O vii and O viii 
lines from sets of observations with the same target pointings. 
They find no universal association between enhanced SWCX 
emission and the closeness of the the line-of-sight to the sub- 
solar region of the magnetosheath. In this paper we do see a 
tendency for XMM-Newton to be clustered around the sub-solar 
region for the SWCX cases, as discussed in Section[3] 

One of the lHenlev & SheltonI (1201 Oh SWCX-enhanced cases 
is part of the data set used in this paper, however this was not 
detected by our method as no discernable variability occuiTed 
during the observation. Cases of SWCX occuiTing in the Earth's 
exosphere, identified by detect ing time-var iable emi ssion, have 
also been observed by Suz aku ( lFuiimoto"et al. 2007: lEzoe et al.l 
I20T0I) . ICarter et all (|20T0^ used a previous observation of a tar- 
get field to constrain the diffuse X-ray emission inherent to that 
look direction to calculate the strength of SWCX emission lines 
associated with a CME passing in the vicinity of the Earth. The 
data set presented in the present paper contains a number of pre- 
viously known cases of SWCX, however, not all of these were 
time variable and therefore were not picked up by the technique 
used here. They did however become tests of the ability of this 
method to identify SWCX-affected observations. One limitation 
of this technique is that SWCX emission occuiTing at an ap- 
proximately steady state as part of any quiescent geocoronal X- 
ray emission will not be identified. In addition, any observed 
quiescent heliospheric SWCX will in general be several times 
stronger than any quiescent geocoronal X-ray emission, due to 
the increased integration lengths involved (Cravens et al. 20011). 
Although a combination of techniques would be ideal, this is not 
possible for the vast majority of single XMM-Newton pointings. 
The greatest advantage of this method is that observations are 
considered on an individual basis. An XMM-Newton user can 
make a judgement as to whether extra caution is required when 
analysing their results for possible SWCX-contamination. 



3.3. Cometary emission 

Comet Hyakutake was the first comet whose X-ray emission, as 
observed by ROSAT and RXTE (Lisse et al. 1996 ), was assigned 
to the SWCX emission process (Cravensi ir997l) . This SWCX 
emission occurs from the interaction of the solar wind with neu- 
tral species that outgas from the comet as it enters the inner so- 
lar system, and the amount of out-gassing is dependent on the 
comet's distance from the Sun. These neutral species are mainly 
water and its dissociation products. The SWCX emission must 
occur in cometary regions where photoionisation and destruction 
of the neutral species can occur. Water and hydroxyl ions have 
short lifetimes when exposed to solar UV photons and therefore 
survive the longest in the coma interior, whereas the dissocia- 
tion products of water (along with CO providing the comet has a 
sufficiently high carbon abundance) can survive further into the 
outer coma regions. A detailed description of X-ray emission 
from comets, prima rily using data from the Chandra obser va- 
tory, can be found in lBodewits et alj (l2007l) : lBodewitsl (l2007l) . 

Several XMM-Newton observations of comets were included 
in the sample in this paper. The x]j and values for the comets 
are given in Table |2] The highest overall values of x^, and 
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Table 2. XMM-Newton observation of comets within the data 
set. 



Revn. 


Obsn. 


Name 


4 




0209 


0103460901 


McNaught-Hartley 


1.5 


1.0 


0369 


0103461101 


C2000 WMl (LINEAR) 


856.7 


299.9 


0719 


0161760101 


Comet 2p (Encke) 


2.0 


1.5 


0720 


0113041301 


C2001 Q4 (Neat) 


1.2 


0.9 


0808 


0164960101 


C2001 Q4 (Neat) 


226.6 


349.77 


1177 


0406950201 


Comet 73p 


13.3 


17.2 



occurred during observations of comets. Example line-band and 
continuum lightcurves from comet C2001 Q4 (Neat) (observa- 
tion 0164960101) are shown in Figure |5] where the line-band 
lightcurve clearly dominates and is highly variable. Also, an 
earlier observation of the same comet exhibited very different 
values of and 7?^. Although the upwind solar wind moni- 
tor Advanced C omposition Explor er (ACE level 2, combined 
instrument data, ^ Stone et aP (Il998h ) proton flux, also shown in 
this graph, is steady and not remarkable in intensity, the comet 
will more likely be sampling solar wind that originates from a 
different location in the solar corona. In Figure |5] we have also 
plotted the to 0*"+ ratio, using data from the ACE SWICS 
instrument. This ratio varies slightly over the length of the obser- 
vation but any fluctuations are not reflected in the XMM-Newton 
line-band lightcurve, further implying that the comet is sampling 
a different solar wind to that seen by ACE. The solar wind is 
a collisionless plasma and so its ion composition remains un- 
changed as it flows away from the Sun. Signatures of S WCX oc- 
curring throughout the solar system can therefore be used to in- 
fer the composition of the solar wind, which varies considerably 
throughout the solar cycle and with solar latitude. As cometary 
orbits are not restricted to the ecliptic plane they are ideal loca- 
tions to study compositional signatures from solar wind originat- 
ing from a variety of solar wind latitudes (Dennerl et al. 1997). 
iBodewits et alj (l2007l) was able to use cometary X-rays to dis- 
tinguish emission resulting from three solar wind types: the cold 
and fast wind, the warm and slow wind and the warm and dis- 
turbed wind. A more complete discussion of cometary X-ray 
emission is beyond the scope of this paper. 



3.4. Planetary emission 

We include within the SWCX set an observation of the planet 
Saturn, taken on 1" October 2002. This observation has been 
comprehensively analysed bv lBranduardi-Ravmont et aH (1201 Ol) 
who attribute the X-ray flux to emission from the planetary disk, 
produced by the scattering of solar X-rays and an additional flu- 
orescent emission line of oxygen at ~0.53 keV origi nating from 
the rin gs. This observation was previously studied bv lNess et alj 
(l2004l) but they did not investigate any low-energy variability. 
We find a high x]i, (2.9) and R.^ (2.3) in our time variability 
test and a distinct step in the line-band lightcurve, indicative of 
a SWCX enhancement, after the same filtering steps as to all 
other data sets have been applied, including source exclusion 
to remove emission from the planet and planetary exosphere. 
The initial source exclusion radius equates to ~ 'i.l Rsatum- 
Spectra from this observation were extracted from event files 
that had been additionally filtered with a larger extraction region 
of ~ \Q.6Rsatum- Even after this additional source extraction 
step the lightcurve production procedure yields metrics of and 
R^ of 3.1 and 2.4 respectively. Two later observations of Saturn 
taken in 2005 were included within our sample but showed no 
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Fig. 5. Lightcurve from comet C2001 Q4 (Neat) that resulted 
in the highest x^fi values (black - SWCX band, red - contin- 
uum band). The time axis is given in hours since the start of 
the XMM-Newton lightcurves. The solar wind proton flux as 
recorded by ACE is given in blue. The O^^ to O^^ ratio is plotted 
using the proton axis and is shown by the blue dashed line. 

evidence of SWCX enhancement. Therefore, we have no rea- 
son to reject this observation from the SWCX set. We were con- 
cerned that the planet's movement through and possibly out of 
the field of view may have caused this effect. However, the dom- 
inant movement is in right ascension with a maximum speed of 
8 arcseconds per hour (Branduardi-Ravmont et al. 2010), result- 
ing in a shift of only 0.78 arcminutes over the course of the 21 ks 
observation. In addition, resolved sources have been removed 
from the field of view as for all other observations. 



4. Spectral analysis 

We continue our study by observing the spectral signatures of 
new SWCX cases identified in this paper and in Paper I. We 
omit those SWCX observations that have been comprehensively 
investigated in the literature yet show no temporal variability in 
our tests. We also omit observations of comets. This set of ob- 
servations which we have used for further study throughout this 
paper is hereafter known as the SWCX set and comprises 103 
observations. 

The Science Analysis System (SAS) software (version 9.0.0; 
http://xmm.esac.esa.int/extemal/xmmjiatajinalysis/) was used 
to produce spectral products and instrument response files for 
all observations of the SWCX set. The SAS accesses instru- 
ment calibration data in so-called cuiTent calibration files (CCFs) 
which are generally updated separately from SAS release ver- 
sions. In this paper we used the latest public CCFs released as of 
February 2010. 

For each exospheric-S WCX case we extracted spectra for the 
EPIC-MOS cameras for the suspected SWCX-affected period 
and for the suspected SWCX-free period. The SWCX-affected 
period was when the enhancement in the line-band lightcurve 
was judged (however not by a formal mathematical argument) 
to have occurred. The enhancement could have occurred at the 
beginning, middle or end of the observation. The remaining time 
periods in the observation made up the SWCX-free period. We 
used events from a circular extraction region, centred on detector 
coordinate positions (DETX, DETY = -50, -180), with an extrac- 
tion radius of 16000 detector units or 13.3 arcminutes. We also 
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Table 3. Principle ion species emission lines used in the model, 
plus any minor emission line energies used (see text). 



Ion 


Energy (keV) 


Minor energies (keV) 


Cv 


0.299 


0.304, 0.308, 0.354, 0.379 


Cvi 


0.367 


0.436, 0.459, 0.471 


Nvi 


0.420 


0.426, 0.431, 0.523 


Nvii 


0.500 


0.593, 0.625, 0.640, 0.650 


Ovii 


0.561 


0.569, 0.574, 0.713, 0.666, 0.698, 0.723 


VIII 


0.653 


0.775,0.817, 0.837, 0.849 


Nex 


1.022 




Mg XI 


1.330 




Si XIV 


2.000 





applied the flag and pattern selection expression '#XMMEA_EM 
&& PATTERN<=12 && FLAG==0'. This pattern selection se- 
lects events within the whole valid X-ray pattern library for the 
EPIC-MOS and the flag selection removes events from or ad- 
jacent to noisy pixels and known bright columns. We produced 
instrumental spectral response files for each period. The instru- 
ment effective area files were calculated assuming the source 
flux is extended, filling the field of view and with no intrinsic 
spatial structure. 



4.1. Spectral modelling 

We knew from the work of ICarter et al.l (1201 Oh that exospheric 
SWCX can occur throughout the entirety of an observation, 
although the line-band lightcurve may show an enhanced and 
a steady-state period. We used the spectra from the apparent 
SWCX-affected period as the source spectra and that from the 
apparent S WCX-free period as the background to produce a dif- 
ference spectrum. The difference spectrum therefore provides a 
lower limit to the SWCX enhancement that has occurred dur- 
ing an observation. Providing the particle-induced background 
is reasonabl y constant over the duratio n of the observation (at 
most +10% i De L uca & Molendil2604 ')) this factor will be ehm- 
inated. An inspection of the difference spectra was made for en- 
ergies above 2.5 keV, to check for the presence of significant 
variable residual soft-proton contamination. Each SWCX-case 
showed a count rate statistically consistent with zero above this 
energy. 

We modelled the resulting difference spectrum for each 
SWCX case with a standard model of emission lines. The spec- 
trum from MOSl and MOS2 were fitted simultaneously, al- 
though a global normalisation parameter for the MOS2 spectrum 
was allowed to vary. The relative line strengths for a particular 
ion species below 1 ke V, for example O vii (which involves seven 
separate transitions including the O vii triplet), were set using the 
velocity dependent cross-sections of laboratory charge exchange 
collisions b etween highly c harged ion and atomic hydrogen, 
as found in iBodewitsI (12007 ). We assumed a solar wind speed 
of 400kms~' for these cross-sections. We also added emission 
lines from Nex at 1.022keV, Mgxi at 1.330keV and Sixiv at 
2 keV. There may be emission from other ion species present in 
the spectra, such as from highly charged iron or aluminium (as 
seen in lCarter et al.l(l2010l) l but we wished to simplify the model 
applied to a general case and the dominant SWCX emission lines 
are found below 1 keV. We fixed the relative normalisations of 
the minor transitions to that of the principal transition for each 
ion species. The principal, dominant transition in the case of C v, 
N VI and O vii is the forbidden line transition. The principal ion 
transitions used in this modelling can be seen in Table |3] 
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Fig. 6. Histogram of total spectrally fitted flux between 0.25 and 
2.5 keV for the SWCX set. 

We used Version 12.5.0 of the XSPECflX-ray spectral fitting 
package to perform this analysis. We fitted the model to each dif- 
ference spectrum by minimising the;i^^ statistic. We calculated 
the modelled flux and 1-sigma errors on the flux between 0.25 
and 2.5 keV for each EPIC-MOS instruments. The keyword giv- 
ing the area collected in the spectral extraction (BACKSCAL) 
was converted into units of steradians and used to convert 
the individual flux values to units of keVcm"^ s"' sr"'. Fluxes 
presented here onwards are for a combined eiTor weighted- 
average EPIC-MOS flux. The error on the flux was calculated 
from the individual flux errors, combined in quadrature. A his- 
togram of the total spectrally fitted flux for each of the SWCX 
set can be seen in Figure |6l The minimum flux we observed 
for a SWCX case was 2.2keVcm"' s"' sr"' (observation id. 
0112490301) and the maximum 50 1 keV cm'^ s"' sr"' (obser- 
vation id. 0085150301. Carter et alj(l2010l) ). 

The relative strengths of the component lines to the SWCX 
spectral model varied considerably within the SWCX set. 
Individual line fluxes were calculated by finding the best fit 
model using all lines, as described above, then setting other line 
normalisations to zero in XSPEC. The flux for the individual ion 
species contribution was calculated in the range 0.25 to 2.5 keV. 

High O^^ to O^^ and magnesium to oxygen ratios are used 
as i ndicators of the presence of CME plasma ( Zhao et al. 2009i 
Ric hardson & Cane 2004; Z urbuchen & Richard son 2006). O vii 
is the dominant SWCX ion-species in the majority of cases. In 
Figure |7] we plot the ratio of the fluxes of the lines Mg xi/O vii 
to Oviii/Ovii (using those oxygen transitions available to us 
within our X-ray spectral band). To look for plasma signatures 
with the highest charge states we only plot those cases where 
the normalisation of the numerator in the ratio is well con- 
strained. Three observations are constrained to have both a ra- 
tio of Mgxi/Ovii > 0.6 and Oviii/Ovii > 1.0. One of these 
(with identifier 0085150301) was the obser vation previously as - 
signed to a passing CME and described in ICarter etalJ (l2010l) . 
These observations are therefore possible candidates for hav- 
ing observed CME plasma with XMM-Newton and are listed 
in Table |4] We quote the lower limit to the ratio in the case 
where the O vii flux is badly constrained, i.e. very weak. In this 
table we also quote the mean value of the O^^ to O^^ ratio dur- 
ing the period of the observation, using values taken from the 
ACE SWICS instrument. This data was only available in two 
of the cases in the table. As ACE is found at Lagrangian point 
LI, we have time shifted the solar wind data to account for the 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 



8 



J. A. Carter et al.: Identifying XMM-A'ewton observations affected by solar wind charge exchange - Part II 



Table 4. SWCX set observations exhibiting the highest 
Mg xi/O VII and O viii/0 vii ratios, or the lower Hmit (95 % con- 
fidence) to this ratio when O vii is badly constrained. We also 
note the ACE SWICS mean value of the to ratio when 
available, with the standard deviation of this ratio given as the 
error. 



Revn 


Obsn 


Mg xi/ 


viii/ 


Mean 






Ovii 


Ovii 


Q7+/06+ 


0342 


0085150301 


2.5±1.9 


8.3±6.4 


0.58+0.54 


0494 


0109120101 


>0.60 


>2.26 


1.53+0.83 



0747 0200730401 >1.09 



>2.49 



travel time to Earth, based on the mean speed of the solar pro- 
tons during the XMM-Newton observation. Expected O^^/O*'^ 
ratios for the slow and fast sol ar wind are 0.27 and 0.03 respec- 
tively (ISchwadron & Cravens! 12000) . Both of the observations 
with ACE SWICS data surpass both the nominal slow and fast 
values by a considerable margin and would suggest that ACE 
detected a CME plasma. Although the identification of CME 
plasma generally involves many more criteria to be satisfied, 
XMM-Newton could provide supplementary spectral evidence 
to studies employing in-situ dedicated solar wind monitors in 
the field of solar system space science. 

To test if any relationship exists between the flux of the 
SWCX lines and increased solar wind flux, we plot in Figure [8] 
the observed flux versus the difference in the mean solar wind 
proton flux (as measured by ACE) between the SWCX-affected 
and SWCX-free periods. We have again time-shifted the ACE 
data to account for the distance between LI and the Earth. 
Although there is considerable scatter amongst these values, 
there is a positive correlation between line flux and solar wind 
proton flux. We include in the plot the linear fit as found by 
the IDL procedure Unfit. Proton flux, for our SWCX set, is a 
good indicator of the presence of SWCX-enhancement, if not 
th e level of this enhancemen t. This is in contrast with the results 
of iHenlev & S helton (2010), whose SWCX cases were consid- 
ered to be due to SWCX occurring within the heliosphere and 
therefore no correlation would be expected between an upstream 
solar wind monitor and any SWCX enhancement. Heliospheric 
SWCX is expected to vary on longer timescales than exospheric 
SWCX and therefore will be harder to identify by the tech- 
nique in this paper. However, at certain times of the year XMM- 
Newton may have a line-of-sight that passes through the helium 
focusing cone, that could potentially produce a variable signal 
in the line-band that may be detectable by t his technique (vari- 
ations over a few hours, jKoutroumpa et al.ll2007l) ). We discuss 
this possibility further in Section 16.11 The flux variations seen 
within our SWCX set are therefore due to local X-ray emission 
in the vicinity of the Earth. 



5. Example cases of SWCX enhancement 

In this section we comment on the three newly identified exo- 
spheric cases from Table [1] Lightcurves for each, over-plotted 
with the solar proton flux as recorded by ACE, are given in 
Figure|9l The solar proton lightcurves have been adjusted for the 
distance between ACE and the Earth by adding a delay based on 
the distance to ACE and the mean solar proton speed, assuming 
a planar wavefront travelling on the Sun-Earth axis. We also plot 
a difference spectrum for each observation by combining data 
from both EPIC-MOS cameras and over-plot the best fit model 
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Fig. 7. Ratio of Mg xi/O vn to O viii/O vii where available for the 
SWCX set. Where appropriate we mark the lower limit (red). 
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Fig. 8. Observed flux versus mean solar wind proton flux. The 
red dotted line indicates a linear fit to the data. 



to the data, as described in Section |4l Individual fluxes for a se- 
lection of prominent lines are given for each case in Table|5] 

- Observation 0150610101 (revolution 0623) 

The line-band lightcurve shows a period of enhanced count 
rate at the beginning of the observation. The ACE solar pro- 
ton flux is raised at the beginning of the lightcurve and re- 
duces as the lightcurve progresses. The difference spectrum 
exhibits emission at O vii, O viii, along with evidence of car- 
bon emission below 0.5 keV. The flux observed between 0.25 
and 2.5 keV was 20.6 keV cm"^ s"' sr"' . 

- Observation 0054540501 (revolution 0339) 

The line-band lightcurve shows an enhancement during the 
latter part of the observation, which is also observed in the 
ACE solar proton flux. The difference spectrum exhibits 
emission in the oxygen band, along with evidence of carbon 
emission below 0.5 keV. The flux observed between 0.25 and 
2.5 keV was 7.9keVcm"^ s"' sr"'. 

- Observation 01 13050401 (revolution 0422) 

The line-band lightcurve shows a period of enhanced count 
rate at the beginning of the observation. A short enhance- 
ment period is seen in the ACE solar proton flux and the 
overall magnitude of this flux is much higher than the other 
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two cases in this section. The flux observed between 0.25 
and 2.5 keV was 25.9keVcm"^ s"' sr"'. 

6. Modelling the expected emissivity 

The expected X-ray emissivity of SWCX emission from the so- 
lar wind interaction with the magnetosheath can be estimated 
from the integrated emission along the line-of-sight for the ob- 
server. We have developed a model, applicable to local interplan- 
etary space, to calculate this emission. We have not attempted 
to include any contribution from further into the heliosphere, as 
to increase the integration length would lead to greater uncer- 
tainty in the underlying parameters of the solar wind on large 
spatial scales. We assume that the solar wind parameters used 
in the model are approximately constant (excluding the magne- 
tosheath region) along the line-of-sight. The emissivity expected 
(ICravensll2000h is given by the expression: 

Px — arjswri«{g) eVcnT^ s^^ (1) 

where a is the efficiency factor dependent on various aspects 
of the charge exchange such as the interaction cross-section and 
the abundances of the solar wind ions, rjsw is the density of the 
solar wind protons, rj,, is the density of the neutral species and 
{g) is their relative velocity. 

For each observation under study in paper II, we wish to test 
whether any relationship exists between the total SWCX flux 
seen and the theoretical integr ated X-ray emiss ion along the line- 
of-sight (based on Equation l2l ICravensI (|2000|) ). 

F=-^r PxdseV cnT^s-^sr^ (2) 
Jo 

We take data describing the conditions in the solar wind 
from ACE (Level2 processed data, merged instrument data us- 
ing hourly averages) at the time of each observation. We needed 
to apply a delay to the signal received, to account for the sepa- 
ration between ACE and the Earth. This delay will be time vari- 
able and will depend on the speed and orientation of the solar 
wind. However, as a first approximation, we have taken the av- 
erage solar proton speed of the data and assumed a planar wave- 
front travelling anti-sunward perpendicular to the GSE-X axis. 
We calculate the delay required for the wavefront to travel from 
ACE to the Earth. 

Throughout this work we assume a geocentric solar-ecliptic 
coordinate system (GSE), where positive X is directed from the 
Earth to the Sun, positive Y opposes planetary motion and posi- 
tive Z is parallel to the direction towards the north ecliptic pole. 

Then for each time bin of an observation: 

- We extract the solar wind proton velocity and temperature 
from the ACE data, and for these parameters we calculate a 
thermal velocity and average speed, using Equations [3] and 

m 

vrh = (3) 

(8> - ^P^ + ^. (4) 

- We estim ate the posi t ion of the magnetopause, based on the 
model of IShue et al. l (fT998l) . To do this we use information 
regarding the strength and direction of the interplanetary 
magnetic field (B, component). We currently assume that the 
magnetopause shape is symmetrical about the GSE-X axis 
and place the magnetopause standoff distance along this axis. 



- Using the magnetopause location as a guide, we approxi- 
mate the position of the bow shock. We base the shape on 
a simple parabola and calculate the bow shock standoff dis- 
tance using the s olar w ind pressure and the relationship in 
iKhan & CowlevI (119991) . The magnetosheath and bowshock 
together define the magnetosheath region. 

- We create an Earth-centric square image for use in subse- 
quent steps. The side length of the image is 200 Rg. This im- 
age is divided into cells, with side length 0.5 Re- The magne- 
topause shape is projected onto this image. We are able to use 
an image rather than a cube due to the assumption made pre- 
viously regarding the symmetry of the magnetosheath shape 
about the GSE-X axis. 

- We find th e neutral density of hy drogen atoms for each cell. 
We use the l0stgaard et al.l(l2003h model for neutral hydrogen 
density profiles arou nd the Earth , but limit this to a minimum 
density of 0.4cm-3 ('Fahi*1971). 

- We determine the line-of-sight of the XMM-Newton point- 
ing through the grid by extracting the relevant information 
from the ODE and converting the positions and target point- 
ing direction to the GSE coordinate system. 

- We fi nd the velocity and density of the solar wind for each 
cell, dSpreiter et al.l (1 1966.) . K.D. Kuntz private communi- 
cation). As the solar wind passes the bowshock and enters 
the magnetopause its density increases (by about a factor of 
four in the subsolar region, as compared to the unperturbed 
value), and the velocity drops to about one tenth. 

- The value of a is dependent on the abundances of the ion 
species contributing to the charge exchange process, along 
with the cross-section and energy of each interaction with 
the neutral donor in the energy band of interest. The neu- 
tral donor is hydrogen in the geocoronal case. The relative 
abundances found in the solar wind vary considerably with 
solar wind state. The composition of the solar wind gener- 
ally follows abundances seen in the photosphere, but can 
vary by up to a factor of about 2 (fast wind) or 4 (slow 
wind) for elements with first ion isation potential (FIP) be- 
low the Lyman-ff limit of 10.2 eV dRichardson & Canel2004 
and references therein). However, we use the slow solar 
wind abundance s for an ion species w ith respect to oxy- 
gen, as listed in ISchwadron & CravensI (l20bQi) . We use an 
oxygen to hydrogen ratio of 1/1780 for solar wind speeds 
of < 650 km s ' or 1/1550 for speeds above this threshold. 
For this modelling we consider contributions to the emis- 
sion from the principal and minor transitions as described in 
Section 14.11 Cross-sections for charge exchange transitions 
are dependent on solar wind speed. We calculate an a map 
with the same dimensions and binning as that of the Earth 
grid and populate this map with values of a depending on the 
speed of the solar wind, unperturbed outside of the magne- 
tosheath or perturbed inside the magnetosheath as described 
above. 

- We multiply the solar wind velocity, solar wind density and 
neutral hydrogen density together for each of the cells in the 
line-of-sight and multiply this value by the efficiency factor 
a for each cell. This is the emissivity of each cell. 

- We sum all cells in the line-of-sight, accounting for the num- 
ber of cells included in the integral, to give the value of the 
emissivity metric, approximating Equation|2] 

There are some known limitations to this model, such as: 

- There are no magnetosheath cusps (increased density or 
modifications to the velocity of the solar wind specific to 
these regions) included in the Spreiter approximation 
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Fig. 9. Example cases of SWCX enhancement, for (top row) observation id. 0150610101, (middle row) observation id. 0054540501 
and (bottom row) observation id. 0113050401. Each left-hand panel shows the line-band (black) and continuum band (red) non- 
mean-adjusted lightcurves along with the solar proton flux (blue, right-hand y-axis). The split between the SWCX-free and SWCX- 
affected period is indicated by the dashed vertical line. In the right-hand panels we show the combined EPIC-MOS diff'erence 
spectrum and the model fitted to the data for each case (solid line). 

Table 5. Most prominent ion line fluxes for example cases. Fluxes are quoted in units of keVcm"^ s ' sr Upper limits (95 % 
confidence) are given for very weak lines. 



Revn Obsn 



Cvi 



Nvi 



Ovii 



O VIII 



0623 0150610101 
0339 0054540501 
0422 0113050401 



4.86±1.10 2.19±1.48 7.73±0.37 2.42+0.61 
5.61±0.85 <3.50 11.73±0.27 4.74±0.41 

9.05+0.41 <2.80 11.60±0.27 4.36+0.46 
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- The neutral hydrogen has been modelled as spherically sym- 
metrical about the Earth. However, there may be density en- 
hancements in regions of the exosphere. 

- The abundances of the solar wind are not constant, but will 
change for example with the phase of the solar cycle or the 
injection of plasma from a CME. 

- The magnetopause and bow shock stand-off distances have 
been assumed to be on the GSE-X axis, this may not be the 
case. 

- The interstellar neutral density may be significantly different 
from that of the approximate limiting density applied in this 
model. 

We also consider an adapted model, whereby the a cal- 
culated is dependent on the relative abundances of O*^ and 
O^^. We use the ratio of these ionisation states, taken from the 
ACE merged, hourly-averaged data sets to re-calculate the abun- 
dance o f O^^ assuming the initial O' '^ and O^^ abundances as 
found in lSchwadron & Cravensl(l2000[) . We then re-calculate the 
value a and the subsequent line-of-sight flux. These results will 
be referred to as the Model-2 results and will be discussed in 
SectionlOl 

6.1. Modelled emission results 

Modelled lightcurves were produced for each SWCX case when 
there were data available from ACE. We split the modelled emis- 
sion based on the time periods used for the creation of the 
spectra in Section 14.11 The resultant flux is the difference be- 
tween the mean modelled flux during the SWCX-affected and 
the SWCX-free periods. A histogram of the modelled fluxes is 
given in Figure [10] (top panel), along with a scatter plot show- 
ing the observed flux versus the modelled resultant flux for each 
exospheric-SWCX case (bottom panel). In general there is a pos- 
itive correlation between the modelled and observed flux. For 
a few cases the modelled flux is negative. This happens when 
the SWCX-affected period, as determined using an enhancement 
seen in the observed line-band lightcurve, occurred in the op- 
posite period to the maximum expected modelled flux. The en- 
hancement in the observed line-band lightcurve occurred suffi- 
ciently far away in time from any peak seen in the modelled 
X-ray flux lightcurve. 

In Figure [TT] (top-row) we show modelled lightcurves for the 
three top new cases of Table \T\ Contributions from the model 
were only taken for the periods where there were counts in the 
XMM-Newton lightcurves (periods not removed during the fil- 
tering process). Example lightcurves of cases where the mod- 
elled flux is negative are given in the first two panels of the sec- 
ond row of Figure [TT] We also plot a modelled lightcurve when 
SWCX was not detected (below the thresholds forx^ and Rj^) in 
the bottom-right panel of the same figure. In this case the line- 
band lightcurve does not vary significantly. The modelled emis- 
sion in this case is small compared to the SWCX cases presented 
in the other examples. 

We wished to test how well the individual modelled flux 
lightcurve tracked that of the line-band lightcurve for each ob- 
servation. We also wanted to determine the most dominant pa- 
rameter in the modelling of the expected emission. To do this 
we applied principal component analysis to the model versus the 
line-band lightcurve and the model versus the solar wind flux. 
We calculate the correlation matrix between a linear fit to the 
relationship between each pair of values, for each exospheric- 
SWCX case. We calculate the correlation rather than the co- 
variance matrix as the scale ranges of the data differ by a large 
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Fig. 10. Top: histogram of the modelled fluxes. Bottom: ob- 
served flux (0.25 to 2.5 keV) versus the modelled flux for the 
SWCX set. A line (dashed, red) of gradient unity has been added 
to the graph to aid the eye. 



amount and so by using the correlation coefficients we standard- 
ise the data. We use the primary eigenvalue of this matrix to 
calculate the percentage contribution along the assumed linear 
relationship between these lightcurves. Histograms of these per- 
centage contributions can be seen in Figure [12] The histogram 
(left) shows that the X-ray lightcurve is generally correlated 
with the modelled lightcurve, as the first principal component 
percentages are high (with a mean of 73.7 %). The histogram 
(right) also shows high first principal component percentages 
(with a mean of 73.6%), which suggests that in the vast ma- 
jority of cases the model is dominated by the incoming solar 
wind flux. The lowest eigenvalue when comparing the modelled 
emission to the line-band lightcurve occurred for the observation 
with the identification number 0101440401. Modelled and line- 
band lightcurves for this case are shown in Figure [T3] We also 
plot the component lightcurves that make up the total modelled 
lightcurve from within the magnetosheath and from beyond the 
bow shock. The contributions from the magnetosheath region in 
this case dominates the modelled lightcurve. XMM-Newton is 
found anti-sunward of Earth during this observation and so the 
line-of-sight of XMM-Newton passed through the flanks of the 
magnetosheath. This region is less well defined in our model due 
to the approximations of the shape of the bowshock boundary 
and the extrapolation of the values used to perturb both the solar 
wind density and velocity in this region. The overall modelled 
emission was very low for this case, compared with the top-row 
cases of Figure nn 
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Fig. 11. Example modelled (blue, in keVcm s sr , left-hand y-axis) and XMM-Newton line-band (black, cts , right-hand 
y-axis) Hghtcurves. Observations with identifiers (top-left) 0150610101, (top-middle) 0054540501, (top-right) 0113050401 show 
the model lightcurve generally following the shape of the XMM-Newton lightcurve. Observations (bottom-left) 0141150101 and 
(bottom-middle) 0150320201 show the modelled lightcurve peak in a different period to the XMM-Newton lightcurve and (bottom- 
right) 0301410601 is an example from an observation without a SWCX enhancement. Five panels show the split between the 
SWCX-affected and SWCX-free periods (vertical dashed line). 
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Fig. 12. Histograms of the first eigenvalues percentage contribution to the total, for the modelled emission versus the XMM-Newton 
line-band lightcurve (left) and versus solar wind flux (right). 



In Figure [14] (top) we plot the modelled resultant flux ver- 
sus the average length of the line-of-sight through the magne- 
tosheath (between the magnetopause and the bow shock). In 
Figure [14] (bottom) we present the observed flux versus the aver- 
age length of the line-of-sight through the magnetosheath. There 
is no discernible general relationship between the modelled or 
the observed flux with the length of line-of-sight through the 
magnetosheath. 

We investigate when the model and observed fluxes are 
discrepant by calculating the fractional difference between the 
observed and modelled fluxes ((observed-modelled)/observed 
flux). In Figure [15] (top panel) we plot this fractional difference 
versus the maximum solar wind flux during each observation, 
along with a histogram of the fractional difference values. The 



mean of these fractional differences was -1-0.17 and the modal bin 
of the histogram was for values between and 1 . A large propor- 
tion (approximately 60 %) of the modelled cases had a fractional 
difference between -1 and 1. The most discrepant cases occurred 
when the solar wind flux was low (compared to the maximum so- 
lar wind flux of these exospheric-SWCX cases). The solar wind 
plasma flow around the Earth's magnetosheath in these cases has 
been badly described by the model. 

The observation with the largest absolute fractional differ- 
ence had identifier 0041750101. This case was similar to cases 
(bottom-left and bottom-right) of Figure [TT] when the mod- 
elled lightcurve peaked in the alternative (SWCX-free) period to 
the enhancement in the observed line-band lightcurve (SWCX- 
affected). 
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Fig. 13. Example modelled Ughtcurve (blue, left-hand y-axis) 
with the XMM-Newton line-band (black, right-hand y-axis), 
for the case where the first eigenvalue percentage contribution 
was the lowest when comparing the modelled flux and XMM- 
Newton lightcurves. The contribution to the modelled lightcurve 
from the magnetosheath (green-dashed) and region past the bow 
shock (plum-dashed) are also shown. The SWCX-aff'ected pe- 
riod was taken between the vertical dashed lines. 
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Fig. 14. Line of sight length through the magnetosheath versus 
the modelled flux for the SWCX set (top) and the observed flux 
for the SWCX set (bottom). A mean error on the observed flux 
bar is given to the right of the bottom plot. 



We also wished to consider whether the fractional diff'er- 
ence was due to some underlying emission with temporal vari- 
ability occurring in near-helio spheric space, in particular to that 
of the helium focusing cone dweller &Meiedfl974 ). We con- 
sider cases within the SWCX set that occur within 10° of the 
cone' s direction (73.9° ecliptic longitude and -5.6° ecliptic lati- 
tude (IWitte et al.l ll996')). As the integration length for the model 
is relatively short compared to the spatial extent of the helium 
focusing cone and size of Earth's orbit, only those observations 
taken when XMM-A^evvfon is within this region are of impor- 
tance. We find 4 cases within this region. These cases are marked 
in red on Figure [15] A statistical analysis, repeatedly drawing 4 
random cases from the SWCX set, indicates that we obtain an 
average fractional difference for the 4 random cases to be greater 
than that of the 4 helium focusing cone cases 28 % of the time. 
We therefore have no evidence to suggest that temporal vari- 
ability originating in the helium focusing cone is a significant 
component of the observed-to-modelled flux discrepancy. 

We also compute the fractional differences between the ob- 
served and modelled flux values for the Model-2 results. These 
are shown in Figure [15] (bottom panel). The peak of the distri- 
bution lies in the same bin as that of Figure [15] (top panel), al- 
though there is a greater variance seen in the differences. We 
conclude that for the SWCX set cases, no benefit has arisen by 
using a compositionally-variable dependent model as opposed 
to the simple model. We continue our discussion based on the 
simple model results only. 

We split the fractional difference values into two sets; for 
cases where this value is < -1.5 or > 1.5 (bad), or any other value 
(good). In Figure [16] we plot histograms of the mid-observation 
position of XMM-Newton (in GSE coordinates, GSE-X, Y and 
Z) for each observation for the good and bad sets. We per- 
formed a Kolmogorov-Smirnov test on the three pairs of good 
and bad sets. The probabilities that the good and bad sets are 
drawn from the same sample distribution were 0.22, 0.002 and 
0.79 for GSE-X, GSE-Y and GSE-Z respectively, indicating 
that for the GSE-Y coordinates, the good and bad sets are sta- 
tistically different. The good set for the GSE-Y positions are 
skewed towards negative values and there are relatively more 
observations in the bad set in the positive direction. We re- 
peat the test using mid-observation position of XMM-Newton 
expressed in Geocentric Solar Magnetospheric (GSM) coordi- 
nates. These differ from the GSE coordinates as the GSM-Y axis 
is perpendicular to the Earth's magnetic dipole (the X-axis is 
unchanged). The Kolmogorov-Smirnov test results were 0.22, 
0.004 and 0.28 for the GSM-X, GSM-Y and GSM-Z respec- 
tively. The Y-coordinate result remains significant. Therefore we 
postulate that the model is better at describing the conditions 
seen by XMM-A^ewfo« when the Y-coordinate is negative. 

The simplifications used in this model to describe the flanks 
of the magnetosheath in terms of shape, solar wind density and 
velocity may mean that the model is less robust in this region. We 
assumed cylindrical symmetry about the GSE-X axis, however, 
the magnetosheath will be non-symmetrical in shape, suffering 
for example magnetosheath erosion along one side of the ma g- 
netopause (along the dusk side, or GSE-Y (lOwen et al ] l2008h l 
a full discussion of which is beyond the scope of this paper. The 
incoming solar wind is expected from the GSE-Y positive direc- 
tion, determined by the flow of the solar wind along the Parker 
Spiral as it emanates from the Sun. It is in this region that we 
expect the greatest differences in shape from the simplified mag- 
netosheath we have used in our modelling steps and it is here that 
we see the largest absolute fractional differences between the ob- 
served and modelled fluxes. It is clear that although the model 
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Fig. 15. Fractional difference between (top) the observed and 
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Newton is found within the helium focusing cone are marked in 
red. 



can estimate the observed flux within a factor of ~ 2 in approx- 
imately 50 % cases, there are still many occurrences when the 
local physical conditions combine so that the simple model does 
not explain the observed flux adequately. 

7. Discussions and conclusions 

We have identified 103 XMM-Newton observations, 3.4 % of the 
sample studied, when temporally variable SWCX emission was 
present in the data. The method of this paper has been able to 
identify cases of temporally variable SWCX from within a large 
sample of XMM-Newton observations. These cases were taken 
from those observations presenting the highest xft ™d val- 
ues. The corresponding occurrence rate within the sample used 
in Paper 1 was ~6.5 %. The data for this paper covered a wider 
range in time compared to Paper 1. The level of detection can 
be attributed to the reduction in solar activity as this time range 
extended into a period towards solar minimum. There will be 
many more XMM-Newton observations affected by SWCX, ei- 
ther occurring within the exosphere or near-interplanetary space, 
such as within the helium focusing cone, or at the heliospheric 
boundary and undetectable here. SWCX occurring within the 
heliosheath will generally vary over longer periods than exo- 
spheric SWCX and so is more suited to detection by observation- 
to-ob servation compariso n (e.g. obs ervations within the stu d- 
ies of iKuntz & SnowdenI (|l008) and Henlev & SheltonI (l2010l) ). 
Enhancements from the helium focusing cone will produce some 
temporal variation but are strongly constrained by viewing ge- 
ometry. The method presented in this paper is only able to iden- 
tify time-variable SWCX which varies over the length of an ob- 
servation and therefore the level of contamination quoted in this 



paper can only provide a lower limit to the occurrence of ex- 
ospheric SWCX as observed by XMM-Newton. When SWCX 
emission is only slowly varying or constant over an exposure it 
will be undetectable by this method. There will also be cases 
which have slipped detection due to a high percentage of the ob- 
servation data being removed by the flare-filtering process, re- 
sulting in short lightcurves that are excluded from our analysis. 
As increased concentrations of solar wind ions in the magne- 
tosheath are expected to mirror increases in the general flux of 
the solar wind, increased levels of SWCX emission are expected 
precisely when the flux of solar proton increases. If the on-board 
radiation monitors of XMM-Newton detect a dangerous environ- 
ment for the satellite, the science instruments are switched into 
a safe mode which invariably leads to the loss of high SWCX 
emission periods being available for detection within our sam- 
ple. Even after soft-proton flare-filtering has been applied to the 
data, considerable proton-contamination may be present. This 
can result in a significant scatter when plotting either the line- 
band or continuum lightcurve, whilst potentially masking a clear 
enhanced period of S WCX-emission during the observation. The 
level of residual soft-proton contamination may mean that the 
observation is completely rejected by an observer If the user 
does indeed proceed to process the data, the limits presented here 
on xf, and may be useful to guide any further analysis as to 
whether extra caution should be taken to account for potentially 
high levels of time-variable SWCX contamination. 

We have shown that exospheric-SWCX occurs preferentially 
on the sunward side of the magnetosheath, when the line-of- 
sight of the XMM-Newton pointing towards its astronomical 
target of interest intersected the area of strongest expected X- 
ray emission of the exosphere. This occurs during the northern 
hemisphere summer months. However, a considerable fraction 
of the SWCX-affected observations had lines-of-sight that in- 
tersected the flanks of the magnetosheath, where the SWCX X- 
ra y emission is expected to be weaker. The example presented 
in ICarter et alJ (12010,) . along with showing the highest flux of 
the SWCX set, is one such case whereby XMM-Newton was not 
pointing in the region of strongest expected X-ray flux. This sug- 
gests that there are considerable deviations from our current un- 
derstanding of either or both the hydrogen neutral density and 
the perturbation of the solar wind in the flank regions of the 
magnetosheath. A dedicated mission observing SWCX emission 
to probe the magnetosheath would answer many questions re- 
garding the distribution of mass and mass transfer i n the magne- 
tosheath, bowshock and near vicinity of the Earth (ICoUier et al.l 
120101). 

For each time-variable exospheric-SWCX case and EPIC- 
MOS instrument, spectra were created for the SWCX-affected 
and the SWCX-free periods. The resulting difference spectrum 
between the two periods became the spectrum used for further 
spectral analysis. We applied, to each difference spectrum, a 
standardised spectral model of 33 Gaussian lines involving 9 
ion species. We set the relative normalisations between lines for 
transitions for one particular species to the ratios of laboratory 
cross-sections measured for a collisional speed of 400 km s"' 
between ions and atomic hydrogen. A combined EPIC-MOS 
flux was calculated between 0.25 and 2.5 keV for each case. 
The SWCX set showed a large spread in spectrally modelled 
observed flux. Although the mean solar proton flux during the 
SWCX-affected period was not a very good indicator of the level 
of observed flux, there was a positive correlation between these 
two parameters. 

The SWCX set showed a range of spectral characteristics, 
with Ovii and Oviii being the dominant lines. Spectral signa- 
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tares obtained from these XMM-Newton observations, such as 
the ratio between magnesium and oxygen ion species, may be 
complementary to data obtained from in-situ solar wind mon- 
itors in classifying solar wind plasma types. We have studied 
the SWCX set with the largest observed flux in a separate paper 
dCarter et al. 2010) which we attributed to a CME passing by 
the Earth. This case also formed part of the subset that exhibited 
the highest Mg xi to O vii and O viii to O vii ratios. CME plasma 
is compositionally different to steady state solar wind plasma. 
Other phenomena, such as co-rotating interacting regions for ex- 
ample, may include high density pulses of plasma but show spec- 
tral signatures close to canonical solar wind plasma conditions. 

We wished to investigate whether the observed spectrally 
modelled flux could be estimated using a simple model, con- 
structed using data describing upwind solar wind conditions, and 
the orbital and target pointing configuration of XMM-Newton at 
the time of each SWCX-affected observation. We used simple 
models of hydrogen densities about the Earth and the perturba- 
tions of the solar wind within the region of the magnetosheath. 
Approximately 60 % of exospheric-SWCX cases showed an ob- 
served to modelled flux fractional difference between - 1 and 1 . 
Negative values of the modelled flux occurred when the model 
predicted an emission pulse in the alternative time period to that 
assigned as the SWCX-affected period. The largest outliers oc- 
curred when the solar wind flux was at its weakest. The model 
was dominated by the solar wind flux. The presence of the mag- 
netosheath made a large contribution to the modelled emission 
in a few cases. The actual line-of-sight length through the mag- 
netosheath did not have any discernable influence on either the 
observed or modelled flux. The model employed a large param- 
eter space and there are various aspects which are expected to 
have a large uncertainty (such as the calculation of the delay 
from ACE to the Earth, if the solar wind plasma front is tilted 
or the distribution of solar wind flow around the magnetosheath, 
especially in the regions far from the subsolar point). Adapting 
the model to account for changes in the solar wind O^^/O''^ ratio 
did not improve the observed to modelled flux fractional differ- 
ence for the SWCX set overall. We have not accounted for any 
anisotropics in the Earth's exosphere in terms of hydrogen den- 
sity. In addition there was some suggestion that those cases when 
XMM-Newton was found at positive GSE-Y (the dusk side) re- 
sulted in the least well-fitting models, where anisotropics in the 
shape of the magnetosheath may be most apparent. 

The authors, who are members of the XMM-Newton EPIC 
BGWG, intend to provide a list detailing the observations af- 
fected by time-variable SWCX in the near future, as part of the 
BGWG web pages and group activities. 
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Table A.l. Table of the SWCX set observations, ranked by^^ (the reduced-;^f^ to the linear fit between the line-band and continuum 
lightcurves). Also Usted for each case are the revolution number (Revn), observation (Obsn) and the MOSl and MOS2 exposure 
identifiers (Expn Ml and Expn M2) and the ratio of the lightcurve variances (Rx)- 



Ind. 


Revn 


Obsn 


Expn Ml 


Expn M2 






Ind. 


Revn 


Obsn 


Expn Ml 


Expn M2 


xi 




1 


0342 


0085150301 


U003 


U003 


27.2 


10.3 


53 


0997 


0206090201 


SOOl 


S002 


2.5 


2.6 


2 


0209 


0093552701 


SOOl 


S002 


23.0 


4.0 


54 


0690 


0134531701 


S002 


S003 


2.4 


2.4 


3 


1014 


0305920601 


SOOl 


S002 


15.0 


30.9 


55 


1023 


0304531801 


SOOl 


S002 


2.4 


2.2 


4 


0690 


0149630301 


SOOl 


S002 


14.1 


21.6 


56 


0476 


0109661201 


SOOl 


S002 


2.4 


2.6 


5 


0623 


0150610101 


U002 


U002 


13.5 


4.8 


57 


1018 


0212480801 


SOOl 


S002 


2.4 


8.8 


6 


0339 


0054540501 


S002 


S003 


13.2 


22.4 


58 


0871 


0206360101 


SOOl 


S002 


2.4 


1.3 


7 


0422 


0113050401 


SOOl 


S002 


12.7 


12.3 


59 


0689 


0149610401 


SOOl 


S006 


2.4 


1.9 


8 


0151 


0094800201 


SOOl 


S002 


12.6 


7.2 


60 


0859 


0203541101 


SOOl 


S002 


2.3 


2.0 


9 


0657 


0141980201 


SOOl 


S002 


12.0 


8.1 


61 


0354 


0049340201 


SOOl 


S002 


2.3 
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